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ABSTRACT 

Advanced  thermal  management  systems  for  internal combustion  engines  can  improve  coolant 
temperature  regulation  and  servomotor  power  consumption  by  6etter  regulating  the 
combustion  process  with  multiple  electro-mechanical  components.  The  traditional  thermostat 
valve,  coolant  pump  and  clutch-driven  radiatorfan  are  upgraded  with  servomotor  actuators. 
Aiken  the  system  components  function  harmoniously,  desired  thermal  conditions  can  6e 
accomplished  in  a  power  efficient  manner.  In  this  paper,  a  comprehensive  control 
architecture  is  proposed  for  transient  temperature  tracking.  An  experimental  system  has  6een 
fabricated  and  assembled  which  features  a  variable  position  smart  thermostat  valve, 
variable  speed  electric  water  pump,  variable  speed  electric  radiator  fan,  engine  blocf  and 
various  sensors.  In  the  configured  system,  the  steam-based  heat  exchanger  emulates  the  heat 
generated  by  the  engine’s  combustion  process.  (Representative  numerical  and  experimental 
results  are  discussed  to  demonstrate  the  functionality  of  the  thermal  management  system  in 
tracking  prescribed  temperature  profiles. 


1.  INTRODUCTION 

Internal  combustion  engine  active  thermal  management  systems  offer  enhanced  coolant  temperature  tracking  control 
during  transient  and  steady-state  operation.  Although  the  conventional  automotive  cooling  system  has  proven 
satisfactory  for  many  decades,  servo-motor  controlled  cooling  components  have  the  potential  to  reduce  the  fuel 
consumption,  parasitic  losses,  and  tailpipe  emissions  (Brace  et  al.,  2001).  Advanced  automotive  cooling  systems 
replace  the  conventional  wax  thermostat  valve  with  a  controllable  position  smart  valve,  and  replace  the  mechanical 
water  pump  and  radiator  fan  with  electric  and/or  hydraulic  driven  actuators  (Choukroun  and  Chanfreau,  2001).  This 
replacement  decouples  the  water  pump  and  radiator  fan  from  the  engine  crankshaft,  and  this  solves  the  problem  of 
having  over/under  cooling,  due  to  the  mechanical  coupling,  and  parasitic  losses  associated  with  running  mechanical 
components  at  high  rotational  speeds  that  increase  exponentially  (Chalgren  and  Barron,  2003).  An  assessment  of 
thermal  management  strategies  for  large  on-highway  trucks  and  high-efficiency  vehicles  were  described  by 
Wambsganss  (1999).  Chanfreau  et  al.  (2001)  studied  the  benefits  of  engine  cooling  with  fuel  economy  and 
emissions  over  the  FTP  drive  cycle  on  a  dual  voltage  42V- 12V  minivan. 

Cho  et  al.  (2004)  investigated  a  controllable  electric  water  pump  in  a  class-3  medium  duty  diesel  engine  trucks.  It 
was  shown  that  the  radiator  size  can  be  reduced  by  replacing  the  mechanical  pump  with  an  electrical  one.  Chalgren 
and  Allen  (2005)  and  Chalgren  and  Traczyk  (2005)  improved  the  temperature  control,  while  decreasing  parasitic 
losses,  by  replacing  the  conventional  cooling  system  of  a  light  duty  diesel  truck  with  an  electric  cooling  system. 

To  create  an  efficient  automotive  thermal  management  system,  the  vehicle’s  cooling  system  behavior  and  transient 
response  must  be  analyzed.  Wagner  et  al.  (2001,  2002,  2003)  pursued  a  lumped  parameter  modeling  approach  and 
presented  multi-node  thermal  models  which  estimate  internal  engine  temperature.  Eberth  et  al.  (2004)  presented  a 
mathematical  model  to  analytically  predict  the  dynamic  behavior  of  a  4.6L  spark  ignition  engine.  To  accompany  the 
mathematical  models,  analytical/empirical  descriptions  were  developed  to  describe  the  smart  cooling  system 
components.  Henry  et  al.  (2001)  developed  a  simulation  model  of  powertrain  cooling  systems  for  ground  vehicles. 
The  model  was  validated  against  test  results  which  featured  basic  system  components  (e.g.,  radiator,  water  pump, 
surge  (return)  tank,  hoses  and  pipes,  and  engine  thermal  load). 

A  multiple  node  lumped  parameter-based  thermal  network  with  a  suite  of  mathematical  models,  describing 
controllable  electromechanical  actuators,  was  introduced  by  Setlur  et  al.  (2005).  The  proposed  simplified  cooling 
system  used  immersion  DC  heaters  to  emulate  the  engine’s  combustion  process  and  control  components,  with 
nonlinear  control  algorithms,  to  regulate  the  temperature.  In  their  experiments,  the  water  pump  and  radiator  fan  were 
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set  to  run  at  constant  speeds,  while  the  smart  valve  was  controlled  to  trade  coolant  temperature  set  point.  Cipollone 
and  Villante  (2004)  tested  three  cooling  control  schemes  ( e.g .,  closed-loop,  model-based,  and  mixed)  and  compared 
them  against  a  traditional  “thermostat-based”  controller.  Page  et  al.  (2005)  conducted  experimental  tests  on  a 
medium-sized  tactical  vehicle  that  was  equipped  with  an  intelligent  thermal  management  system.  The  authors 
investigated  improvements  in  the  engine’s  peak  fuel  consumption  and  thermal  operating  conditions.  Finally, 
Redfield  et  al.  (2006)  operated  a  class  8  tractor  at  highway  speeds  to  study  potential  energy  saving  and  demonstrate 
engine  cooling  to  with  ±3°C  of  set  point  value. 

In  this  paper,  a  nonlinear  control  strategy  is  presented  to  actively  regulate  the  coolant  temperature  in  internal 
combustion  engines.  An  advanced  thermal  management  system  was  implemented  on  a  laboratory  test  bench  that 
features  a  smart  thermostat  valve,  variable  speed  electric  water  pump  and  fan,  radiator,  engine  block,  and  a  steam- 
based  heat  exchanger  to  emulate  the  combustion  heating  process.  The  proposed  control  strategies  have  been  verified 
by  simulation  and  validated  by  experimental  testing.  In  Section  2,  the  cooling  system  model  is  presented  to  describe 
the  thermal  system  dynamics.  Nonlinear  tracking  control  strategies  are  introduced  in  Section  3.  Section  4  presents 
the  experimental  steam  test  bench  while  Section  5  introduces  numerical  results.  Experimental  results  are  introduced 
in  Section  6.  The  conclusion  is  contained  in  Section  7.  The  Appendices  present  a  Lyapunov-based  stability  analysis, 
which  is  needed  for  the  controller’s  design,  as  well  as  the  Nomenclature  List. 

2.  COOLING  SYSTEM  MODEL 

A  reduced  order  two-node  lumped  parameter  thermal  model  is  used,  as  shown  in  Fig  1,  to  represent  the  transient 
response  of  the  advanced  cooling  system.  The  system  components  are  a  4.6L  engine  block,  a  radiator  block  with  a 
variable  speed  electric  fan,  a  variable  position  three-way  smart  valve,  a  variable  speed  electric  water  pump,  and  a 
steam-based  heat  exchanger  to  emulate  the  heat  generated  by  the  combustion  process. 


Fig  1  Advanced  cooling  system  which  features  a  smart  valve,  variable  speed  pump,  variable  speed  fan,  engine  block, 
radiator,  and  sensors  (temperature,  mass  flow  rate,  and  power) 

This  simplified  model  is  used  to  reduce  the  computational  burden  for  in-vehicle  implementation.  The  dynamics  of 
the  engine  and  radiator  blocks  at  the  selected  nodes,  as  shown  in  Fig  1,  can  be  described  by  the  following 
differential  equations 


CeT=Qln-Cpcmr(T-T)  (1) 

CrTr  =  Cpcmr  ( Te  -Tr)  —  Q0+s Cpam f  (Te  -Tx).  (2) 

Because  a  three-way  valve  is  used  in  the  system,  the  variable  mr  (t)  is  defined  as  mr  =  Hmc ,  where  the  variable 
H(t)  satisfies  the  condition  0  <  H(t)  <  1 .  Note  that  H(t )  =  1(0)  corresponds  to  a  Hilly  closed  (open)  valve  position 
and  coolant  flow  through  the  bypass  (radiator)  loop.  To  facilitate  the  controller  design  process,  four  assumptions  are 
imposed. 

Assumption  1:  The  system  parameters;  Cpc ,  Cpa ,  C, ,  Cr  and  £  are  assumed  to  be  constants  and  fully  known. 
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Assumption  2:  It  is  assumed  that  the  engine  block  and  radiator  block  temperatures  satisfy 
Te(t)  -  Tfit)  >  s2,\/t  >  0  where  s2  e  SJ?+  is  a  constant.  Further  we  assume  that  7^(0)  >  Tj.jO)  to  facilitate  the 
boundedness  of  signal  argument. 

Assumption  3:  It  is  assumed  that  the  surrounding  ambient  temperature  T.,  ( t )  is  constant  and  satisfies 
Tr(t) -Tm(t)  >  £l,  Vt  >0  where  £l  e  iR'  is  a  constant. 

Assumption  4:  From  the  problem  physics,  the  signals  Qin(t)  and  Q0{t)  always  remain  positive  in  (1),  and  (2) 

3.  THERM  AT,  CONTROL  SYSTEM  DESIGNS 

In  this  paper,  two  Lyapunov-based  nonlinear  control  algorithms  are  presented  to  maintain  a  desired  thermal 
condition  for  the  engine  block.  The  objective  of  the  proposed  control  strategies  is  to  get  a  precise  engine  temperature 
tracking  while  compensating  for  the  system  uncertainties  (e.g.  combustion  process  heat  input  and  external  ram  air 
disturbance)  by  harmoniously  controlling  the  system  actuators. 

The  control  objective  of  the  nonlinear  control  strategies  is  to  ensure  that  the  actual  temperatures  of  the  engine  block 
Te(t)  and  the  radiator  block  Tfit)  track  the  desired  trajectories  Ted(t)  and  Tnfit) ,  respectively,  in  the  following 
sense 

Tfit)  -+Ted(t)  as  t  — »  oo  ,  pfit)^Trd(t)  as  ?  — »  oo  (3) 

while  compensating  for  the  system  uncertainties  Qin(t)  and  Qa(t) . 

Remark  1:  It  is  assumed  that  the  desired  engine  and  radiator  temperature  profiles  are  selected  such  that  they  are 
always  bounded  and  their  first  three  time  derivatives  remain  bounded  at  all  times  ( i.e ., 
TeA0,tA0,tA0Jed(t)Jld(t),Tld(t),frd(tl  TrAt)  e  Lm  .  Further  more,  TeAt)  »  Trd(t)  »  Tm  at  all  times. 

3.1  Adaptive  Control  Strategy 

To  facilitate  the  controller’s  development  and  quantify  the  temperature  tracking  control  objective,  the  following 
signals  are  defined 


a  t  n  =  T  -T 

e  1ed  1e  ’  7  rd  1r  • 

Assumption  5:  It  is  assumed  that  the  signals  Qin  (t  )  and  Qa  (/)  are  constants  at  all  time. 


(4) 


3.1.1  Closed-Loop  Error  System  Development  and  Controller  Formulation 

The  open-loop  error  system  can  be  analyzed  by  taking  the  first  time  derivative  of  both  expressions  in  (4) 
and  then  multiplying  both  sides  of  the  resulting  equations  by  Ce  and  Cr  for  the  engine  and  radiator  dynamics, 
respectively.  Thus,  the  system  dynamics  described  in  (1)  and  (2)  can  be  substituted  and  then  reformatted  to  realize 

Cee  =  CeTed  -0e  +  u e ,  Crrj  =  CrTrd  +  6r  -  ur  (5) 

where  9e{f)  =  Qjt),  6r(t)=Q0(t),  ue  (t)  =  Cpcmr  (Te  -  Tr ) ,  and  u\t)  =  Cpcm\Te-Tr)-£Cpamf(Te-TA). 

Remark  2:  The  control  inputs  mr  (f)  and  rhAt)  are  uni-polar.  Hence,  commutation  strategies  are  designed  to 
implement  the  bi-polar  inputs  ue(t )  and  ur(t )  as 

.  A  ^[i  +  sgn^J]  .  A  F[l  +  sgn(F)] 

m= — - — t - f1-,  mf= — - t - 4-  (6) 

2Cpc{Te-Tr)  f  laCjT-A) 

where  F(t)  =  Cpcmr  (Te -Tr^-ur .  The  control  input,  mf(t)  is  obtained  from  (6)  after  mr(t )  is  computed.  From 
these  definitions,  it  is  clear  that  if  ue  ( t),ur  (t)  g  Lfis/t  >  0 ,  then  mr  ( t),mf  (t)  e  Lfs/t  >  0  .  From  the  calculation  of 
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the  radiator  mass  flow  rate  »/,.(?)  and  using  a  three-way  valve  in  the  system,  the  water  pump  speed  can  be 
determined  for  a  given  valve  position  or  the  valve  position  can  be  determined  for  a  given  water  pump  speed. 


To  satisfy  the  control  objectives  stated  earlier,  continuous  control  laws  ue(t)  and  «,  (/  )can  be  designed  as  follows 

we  =  &e  -  CeTed  -  Kce  ,  Ur  =  8r  +  CrTrd  +  K,n  (7) 


where  the  estimates 6e{t)  and  9r{t)  are  to  compensate  for  the  system  constants  uncertainties  Qin ( t )  and Q0{t) ,  and 
designed  as  follows 


6e  =  ~ae  f  e(f)  dT  >  dr=<Xrljl{*)dT 


(8) 


3.1.2  Sta  bility  Analysis 

Theorem  1:  The  controller  given  in  (7)  ensures  asymptotic  engine  and  radiator  temperatures  tracking  (i.e. 
e{t\rj(t)  — >  0  as  t  — >  oo  )  and  all  closed-loop  signals  are  bounded. 

Proof:  See  Appendix  A  for  the  complete  Lyapunov-based  stability  analysis. 

3.2  Robust  Control  Strategy 

To  facilitate  the  controllers’  development  and  quantify  the  temperature  tracking  control  objective,  a  fdtered 
tracking  error  is  defined  as  follows 


se  =e  +  aee  ,  sr=p  +  arp 


(9) 


where  eft)  and  ij(t)  were  defined  in  (4). 


Assumption  6:  It  is  assumed  that  the  signals  (Qln  (t)  and  Q„(t)  vary  with  time  and  their  first  two  time 
derivatives  remain  bounded  at  all  time,  such  that  Qm  (t),  Q,„  (t),  Qm  (t),  QJt),  Q0  (t),  Q0  (t)  e  Lm  . 


3.2.1  Closed-Loop  Error  System  Development  and  Controller  Formulation 

The  open-loop  error  system  can  be  analyzed  by  taking  the  first  time  derivative  of  both  expressions  in  (9)  and  then 
multiplying  both  sides  of  the  resulting  equations  by  Ce  and  Cr  for  the  engine  and  radiator  dynamics,  respectively. 
Thus,  the  system  dynamics  described  in  (1)  and  (2)  can  be  substituted  and  then  reformatted  to  realize 


CJe  =  CeTed  ~  Qin  +  . 


C  s  =  C  T  .  +  Q  —  u  +  C  a  it 

r  r  r  rd  r  r  r  ' 


(10) 


where  (4)  and  it  is  first  time  derivative  were  utilized  provided  the  signals  ue(t )  and  u,.{t)  are  defined  as 
«,  (0  =  Cpcmr  (Te-Tr)  and  ur  (/)  =  Cpcmr  (Te-Tr)-e  Cpamf  (T -Tr). 

Remark  3:  The  control  inputs  mr  (t)  and  mf(t )  are  uni-polar.  Hence,  the  commutation  strategies,  designed  in 
(6),  are  used  to  implement  the  bi-polar  inputs  ue  (t)  and  ur  (t) . 


To  facilitate  the  subsequent  analysis,  the  expressions  in  (10)  are  rewritten  as  follows 

Cese  =Ne+  Ned  +ue  -  e  ,  Crsr  =Nr+  Nrd  -  ur  -  p 

where  the  auxiliary  signals  Ne(fe,Te,t)  and  Nr(fr,Tr,t )  are  defined  as  follows 


(11) 


N=N-N„ 


N=N-  N„ 


(12) 
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where  Ne{fe,Te,t)  and  Nr (fr , Tr , t)  are  defined  as  follows 

K  =  cefed  -  Qin  +  Ctaee  +  e,  Nr=  CrTrd  +  Q0  +  Crarrj  +  p 
and  both Ned(t)  and  Nrd(t)  are  defined  as  follows 


Ned  ~Ne  — CeTed  Qin, 


Nrd  =  Nr  j >Trd  jr=irJ  ~  CrTrd  +Qo  ■ 


Based  on  (12)  through  (14),  the  control  laws  introduced  in  (11)  are  designed  as  follows 

ue  =  ~(Ke  +  «J[e  -  eJ-J'  [a,{Ke  +  ae)e(r)  +  pe  sgn(e(r))]  d r 
ur  ={Kr+  a, )[/;  -  77(,]+  ['  [a,  ( A',  +  ar)Tj(T)  +  pr  sgn(/;(r))]  dr 


(13) 

(14) 

(15) 

(16) 


where  eo,i)o  are  e{t\ri{t) ,  respectively,  computed  at  the  initial  time  to .  The  terms  ea  and  qo  in  (15)  and  (16), 
respectively,  are  included  so  that  ue(ta)  and  u,,(t0)  equal  to  zero,  where  the  “sgn”  terms  compensate  for  the 
unknown  quantities  in  (13).  The  time  derivatives  of  (15)  and  (16)  are  given  by  the  following  expressions 


ue  =  ~{Ke  +  ae)se  -  pc  sgn(e) ,  ur  =(Kr+  ar ) .s;  +  pr  sgn(/7) . 

After  substituting  (17)  into  (11),  the  closed-loop  error  system  can  be  obtained  as  follows 

Cese  =Ne+  Ned  -  ( Kc  +ae)se-pe  sgn(e)  -  e 
C,.sr  =  Nr  +  Nrd  -  (Kr  +  ar  -  p,.  sgn(/7)  -  rj . 


(17) 


(18) 

(19) 


Remark  4:  Based  on  Remark  1,  Assumption  6,  and  the  expressions  in  (14),  (?)|, | (?)|, (? )|  and  |Arrf(t)| 

can  be  upper  bounded  by  known  positive  constants  as  follows 


\Ned\<Ce  1, 


Ned  ^2, 


Nrd  <Cr2- 


(20) 


3.2.2  Stability  Analysis 

Theorem  2:  The  controller  given  in  (15)  and  (16)  ensures  asymptotic  engine  and  radiator  temperatures  tracking 
(i.e.  e(t\se(t\rt{t\sr(t)^>  0  as  t  — >  00  )  and  all  closed-loop  signals  are  bounded  provided  the  control  gains  pe  and 
pr  are  selected  to  satisfy  the  following  sufficient  conditions 

Pe>£el+—  Ce2,  Pr  >  Crl+ ~Cr2  (21) 

ar 

where  ,  ^e2 , and  If rl  are  given  in  (20),  and  Ke  and  Kr  are  selected  sufficiently  large. 

Proof:  See  Appendix  B  for  the  complete  Lyapunov-based  stability  analysis. 

4.  STEAM  TEST  BENCH 

An  experimental  test  bench  (refer  to  Figure  2)  has  been  fabricated  to  demonstrate  the  proposed  advanced  thermal 
management  system  controller  design.  The  assembled  test  bench  offers  a  flexible,  rapid,  repeatable,  and  safe  testing 
environment.  Clemson  University  facilities  generated  steam  is  utilized  to  rapidly  heat  the  coolant  circulating  within 
the  cooling  system  via  a  two-pass  shell  and  tube  heat  exchanger.  The  heated  coolant  is  then  routed  through  a  6.0L 
diesel  engine  block  to  emulate  the  combustion  process  heat.  From  the  engine  block,  the  coolant  flows  to  a  three-way 
smart  valve  and  then  either  through  the  bypass  or  radiator  to  the  water  pump  to  close  the  loop.  The  thermal  response 
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of  the  engine  block  to  the  adjustable,  externally  applied  heat  source  emulates  the  heat  transfer  process  between  the 
combustion  gases,  cylinder  wall,  and  water  jacket  in  an  actual  operating  engine.  As  shown  in  Figure  1,  the  system 
sensors  include  three  type-J  thermocouples  (e.g.,  T[  =  engine  temperature,  T2  =  radiator  temperature,  and  T3  = 
ambient  temperature),  two  mass  flow  meters  (e.g.,  Mi  =  coolant  mass  flow  meter,  and  M2  =  air  mass  flow  meter), 
and  electric  voltage  and  current  measurements  (e.g.,  Pi  =  valve  power  consumed,  P2  =  pump  power  consumed,  and 
P3  =  fan  power  consumed). 


Fig  2  Experimental  thermal  test  bench  that  features  a  6.0L  diesel  engine  block,  three-way  smart  valve,  electric 
water  pump,  electric  radiator  fan,  radiator,  and  steam-based  heat  exchanger 

The  steam  bench  can  provide  up  to  55  kW  of  energy.  High  pressure  saturated  steam  (412  kPa)  is  routed  from  the 
campus  facilities  plant  to  the  steam  test  bench,  where  a  pressure  regulator  reduces  the  steam  pressure  to  172  kPa 
before  it  enters  the  low  pressure  filter.  The  low  pressure  saturated  steam  is  then  routed  to  the  double  pass  steam  heat 
exchanger  to  heat  the  system’s  coolant.  The  amount  of  energy  transferred  to  the  system  is  controlled  by  the  main 
valve  mounted  on  the  heat  exchanger.  The  mass  flow  rate  of  condensate  is  proportional  to  the  energy  transfer  to  the 
circulating  coolant.  Condensed  steam  may  be  collected  and  measured  to  calculate  the  rate  of  energy  transfer.  From 
steam  tables,  the  enthalpy  of  condensation  can  be  acquired.  To  facilitate  the  analysis,  pure  saturated  steam  and 
condensate  at  approximately  7’=100°C  determines  the  enthalpy  of  condensation.  Baseline  testing  was  performed  to 
determine  the  average  energy  transferred  to  the  coolant  at  various  steam  control  valve  positions.  The  coolant 
temperatures  were  initialized  at  Te  =  67°C  before  measuring  the  condensate.  Each  test  was  executed  for  different 
time  periods. 

5.  SIMULATION  RESULTS 
5.1  Adaptive  Control 

Test  1:  Constant  Heat  Input  and  External  Disturbance 

A  numerical  simulation  for  the  adaptive  control  strategy,  introduced  in  Section  3.1,  has  been  performed  on  the 
system  dynamics  described  in  (1)  and  (2)  to  demonstrate  the  performance  of  the  proposed  controller  given  in  (7)  and 
(8).  The  simulated  thermal  system’s  parameters  were  set  to  the  following:  Ce  =  Cr  =  63250  [J/°K], 

Cpc  =4186  [J/kg.°K],  Cpa  =  1000  [J/kg.°K],  £■  =  0.6  .  T„(t)  =  300  [°K],  rhc(t)=  0.8214  [kg/sec],  Qj$)  =  30  [kW] 
and  Q0(t)=  7  [kW].  The  initial  simulation  conditions  were  set  to  be  as  follows:  Te(o)=  360  [°K]  and 
7j(o)  =  347  [°K],  The  desired  temperatures  were  selected  to  be  Ted(t)  =  350  +  sin(0.05t)  [°K]  and  Trd{t)  =  340  [°K]. 
The  controller  gains  were  set  to  be  as  follows:  Ke  =  1000  ,  Kr  =  10  and  ae  =  ar  =  0.1 .  The  engine  and  radiator 
temperatures  response  is  presented  in  Fig  3,  while  the  normalized  valve  position  and  fan  response  is  shown  in  Fig  4. 
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Fig  3  Engine  &  radiator  temperatures  response  Fig  4  Normalized  valve  position  &  fan  response 


Test  2:  Variable  Heat  Input  and  External  Disturbance 

A  numerical  simulation  for  the  adaptive  control  strategy,  introduced  in  Section  3.1,  has  been  performed  on  the 
system  dynamics  described  in  (1)  and  (2)  to  demonstrate  the  performance  of  the  proposed  controller  given  in  (7)  and 
(8).  The  simulated  thermal  system’s  parameters  were  set  as  in  Test  1  except  for  Qh,(t)  and  Q0{t),  where 
(2o(t)=10  +  3sin(0.05t)  [kW]  and  Q,„{t)  as  shown  in  Fig  5.  The  heat  input  is  designed  as  shown  in  Fig  5  to  verify 
the  reliability  of  the  control  design  and  performance  under  different  levels  of  disturbance.  The  initial  simulation 
conditions  and  desired  temperatures  were  also  selected  as  in  Test  1.  The  controller  gains  Kr  and  a,,  were  reset  to  be 
as  follows:  Kr  =  500  and  a,.  =  0.01 .  The  engine  and  radiator  temperatures  response  is  presented  in  Fig  6,  while  the 
normalized  valve  position  and  fan  response  is  shown  in  Fig  7. 
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Fig  6  Engine  and  radiator 
temperatures  response 


Fig  7  Normalized  valve  position 
and  fan  response 


When  the  heat  input  is  set  to  be  high  (40  or  35  [kW]),  the  valve  routs  the  coolant  only  through  the  radiator  loop  and 
the  fan  runs  at  its  maximum  speed  as  shown  in  Fig  7.  The  valve  and  fan  efforts,  in  fact,  are  not  enough  to  cool  down 
the  engine  temperature  as  shown  in  Fig  8  because  of  the  extra  heat  added  to  the  system  that  can  not  be  handled. 


5.2  Robust  Control 

Test  3:  Constant  Heat  Input  and  External  Disturbance 

A  numerical  simulation  for  the  robust  control  strategy,  introduced  in  Section  3.2,  has  been  performed  on  the  system 
dynamics  described  in  (1)  and  (2)  to  demonstrate  the  performance  of  the  proposed  controller  given  in  (15)  and  (16). 
The  simulated  thermal  system’s  parameters,  initial  simulation  conditions,  and  desired  temperatures  were  set  as  in 
Test  1.  The  controller  gains  were  set  to  be  as  follows:  Ke  =  1000  ,  Kr  =  100  ,  ae  =  a,.  =  0.1  and  pe  =  pr  =  0.5  .  The 
engine  and  radiator  temperatures  response  is  presented  in  Fig  8,  while  the  normalized  valve  position  and  fan 
response  is  shown  in  Fig  9. 
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Fig  8  Engine  &  radiator  temperatures  response 
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Fig  9  Normalized  valve  position  &  fan  response 


Test  4:  Variable  Heat  Input  and  External  Disturbance 

A  numerical  simulation  for  the  robust  control  strategy,  introduced  in  Section  3.2,  has  been  performed  on  the  system 
dynamics  described  in  (1)  and  (2)  to  demonstrate  the  performance  of  the  proposed  controller  given  in  (15)  and  (16). 
The  simulated  thermal  system’s  parameters,  initial  simulation  conditions,  and  desired  temperatures  were  set  as  in 
Test  1.  The  controller  gains  were  reset  to  be  as  follows:  XE=1000,  Kr  =  500  ,  =  5,  ar  =0.1  and 

pe  =  pt.  =  0.01 .  The  engine  and  radiator  temperatures  response  is  presented  in  Fig  10,  while  the  normalized  valve 
position  and  fan  response  is  shown  in  Fig  11.  From  Fig  10,  it  is  clear  that  the  robust  control  can  not  handle  the  extra 
heat  added  to  the  system  as  in  the  adaptive  control. 
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Fig  10  Engine  and  radiator  temperatures  response 


Fig  11  Normalized  valve  position  and  fan  response 


Table  1  shows  a  summery  for  the  simulation  results.  From  Table  1,  it  is  clear  that  the  steady  state  of  the  engine 
temperature  tracking  error  is  about  the  same  for  all  the  proposed  control  strategies.  It  is  also  clear  that  the 
performance  of  both  control  strategies  are  quite  the  same  in  terms  of  power  measure  and  handling  the  heat 
variations.  The  power  measure  is  found  as  follows 


P,yS=[WM+m)iT<l\dT 


(22) 


where  the  power  consumed  by  the  system  actuators  mainly  depends  on  how  much  coolant  and  air  mass  flow  rates 
are  needed.  On  the  other  hand,  the  power  consumed  by  the  valve  is  considered  very  small  and  neglected. 


Test  # 

Description 

Steady  State  Error  [°K] 

Power  Measure 

1 

Adaptive  control  ( Qin  /  Qn  are  constants) 

[-0.03,  0] 

2952 

2 

Adaptive  control  ( Qin  /  Qa  are  variables) 

[-0.03,  0] 

2529 

3 

Robust  control  ( Qin  /  Qa  are  constants) 

[-0.03,  0.015] 

2939 

4 

Robust  control  ( Qm  /  Qa  are  variables) 

[-0.03,  0.005] 

2528 

Table  1  Simulation  results  summery 


6.  EXPERIMENTAL  RESULTS 


7.  CONCLUSION 


Advanced  automotive  thermal  management  system  can  have  a  positive  impact  on  gasoline  and  diesel  engine  cooling 
systems.  In  this  paper,  a  suit  of  servo-motor  based-cooling  system  components  have  been  assembled  and  controlled 
using  a  Lyapunov-based  nonlinear  control  technique.  The  control  algorithm  has  been  investigated  using  both 
simulation  and  experimental  tests.  Two  detailed  and  two  supplemental  controllers  were  applied  to  regulate  the 
engine  temperature.  In  each  instance,  the  controllers  successfully  maintained  the  engine  block  to  setpoint 
temperatures  with  small  error  percentages.  It  has  also  been  shown  that  the  power  consumed  by  the  system  actuators 
can  be  reduced.  Overall,  the  findings  demonstrated  that  setpoint  temperatures  can  be  maintained  satisfactory  while 
minimizing  power  consumption  which  ultimately  impacts  fuel  economy. 
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APPENDIX  A 
Proof  of  Theorem  1 

Let  V ( z ,  t)  g  SR  denote  the  following  non-negative  function 


V  =  —  C  e2  +  —92  +-C \rf  +  —  9,2 
2  ‘  2a„  2  2a 


where  z(t )  g  SJ12  is  defined  as  follows 


z  =  e  >l\‘  ■ 


Note  that  (22)  is  bounded  as  follows 


^||z(?)f<F(z,t)<T2||z(t)||2 


(22) 


(23) 


(24) 


where  are  positive  constants.  After  taking  the  time  derivative  of  (22),  the  following  expression  can  be 

obtained 


V  =  e[de-  Kee)  +  —  9ede  +  /;(-  9r  -  Krn)+  —  0r0r 
ae  ar 

where  (5)  and  (7)  were  utilized,  and  the  auxiliary  signals  6r (t )  and  0r(t)  are  defined  as  follows 


0^0-0,  0^0-0 

e  e  e  5  r  r  r 


The  expression  in  (25)  can  be  written  as  follows 

V  =  e9e  - K /'  +-0e{- aee) - >{0r  ~ K/j2  +—9r (a, :rj) 


cc. 


(25) 


(26) 


(27) 


where  the  expressions  in  (8)  and  their  first  time  derivatives  were  utilized.  By  utilizing  (23),  V(t)  of  (27)  can  be 
rewritten  and  upper  bounded  as  follows 


V<-X ,  z 


(28) 
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where  A3  =  min {Ke,Kr}.  From  (22),  (24)  and  (28),  it  is  clear  that  V(z,t)  e  Lx  ;  hence,  e(t),ij(t),z(t )  e  .  From  (4) 
and  Remark  1,  it  is  clear  that  Te(t),Tr(t)  e  L ^  .  From  (7),  (8),  and  Remark  1,  it  is  clear  that  ue(t),ur(t )  e  L w  .  From 
(5)  and  the  previous  bounding  statement,  it  is  clear  that  e(t),r/(t)  e  Lx  ;  hence,  Te(t),Tr(t)  e  ,  where  (4)  was 
utilized.  From  (6),  and  also  the  previous  bounding  statement,  it  is  clear  that  mr(t),mf(t)  e  Lm  ;  thus,  it  is  clear  that 
mc(t),H(t)  e  Lm  .  A  direct  application  of  Theorem  8.4  in  (FI.  Khalil,  2002)  can  be  used  to  prove  that  z(t )  — >  0  as 
t  — >  oo  ;  thus,  e(t),?)(t)  — >  0  as  t  —>  oo  ,  where  (23)  was  utilized.  Based  on  the  definition  of  e{t )  and  i](t)  in  (4),  it  is 
clear  that  if  e(t),?j{t)^  0  as  t  — >  oo  ,  then  Te(t)  Telj(t),Tr(t)  Trcl(t)  as  t —>  oo  .  Since  (24)  and  (28)  hold 
globally,  from  Theorem  8.4  in  (FI.  Khalil,  2002),  it  is  clear  that  this  proof  provides  a  global  asymptotic  stability 
(GAS)  result. 


APPENDIX  B 
Proof  of  Theorem  2 

Lemma  1:  Let  the  auxiliary  functions  Le(t),Lr(t)  £  5?  be  defined  as  follows 

Le=se(Ned-pesgn(e)),  Lr=sr(Nrd-prsgn(ji))  (29) 

where  if  pc,pr ,  introduced  in  (15)  and  (16)  respectively,  are  selected  to  satisfy  the  conditions  in  (21),  then 


lLXr)dr<Ceb,  \‘Lr{r)dT<Crb 


(30) 


where  the  positive  constants  Ceb,Crb  G  -h  afe  defined  as  follows 


Pe\e{h)\-e{to)Ned(t0), 


Crb=pr  \ri{t0)\-ri{to)N,.d{t0)  . 


(31) 


Proof:  After  substituting  the  first  expression  in  (9)  into  the  first  expression  in  (29),  and  then  integrating,  the 
following  expression  can  be  obtained 


{  4 (r)dr  =  ae  j'  e{r)[Ned  -  pe  sgn(e(r))]</ r  +  \‘J^-NeddT -  pe  j'  ^^sgn(e(r )) d r  . 


(32) 


After  evaluating  the  second  integral  on  the  right-hand  side  of  (32)  by  parts  and  evaluating  the  third  integral,  the 
following  expression  is  obtained 


[  LXr)dT=ae^  e{T)[Ned  -  pe  sgn(e(r))]t/ r  -  J'  e(r)  dr+  e(r)Ned  (t)|^  -  p\e{r^  .  (33) 


The  expression  in  (33)  may  be  upper  bounded  as  follows 


f  4(r)rfr  =  f  \e{r) 


aeNed  + 


dNed{r)\ 


dr 


-a  p 

ere 


dr  +  \e(t)(\Ned {t)~  pe)+ pe \e(ta )| -  e(t0 )Ned (tj.  (34) 


From  (34),  it  is  easy  to  see  that  if  is  chosen  according  to  (31),  then  the  first  inequality  in  (30)  is  hold.  The  second 
inequality  in  (30)  can  be  obtained  by  integrating  the  second  expression  in  (29)  as  follows 


+ 


dNjj)  I 


dr 


■  arP, 


dr  +  -  pr )+  pr  |  T](to  )|  -  ;7(to  )Nrd  (ta ) .  (35) 
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Proof  of  Theorem  2:  Let  the  following  functions  Pe(t), Pr(t) &  i R  be  defined  as  follows 

Pe=Ceb-\'tLe(T)dT>  0,  Pr^^-J'4(r)rfr^0  (36) 

where  Le(t),Lr{t),£eb,£rb  were  defined  in  (29),  and  (31).  The  results  from  Lemma  1  can  be  used  to  show  that 
Pe(t),  P,{t)  are  non-negative  functions.  Let  V(y,t )  e  51  denote  the  following  non-negative  function 


V  =  -e2+-  C  s  2  +Pe+-if+-  C  s  2  +  Pr 
2  2  2  2 


where  y(t )  e  516  is  defined  as  follows 


(37) 


y  =  [z  ^  (38) 

where  the  composite  vector  z(t )  e  514  is  defined  as  follows 

Z=[Ze  ZrJ=[e  ^  *1  ■Sr]7’  (39) 

where  ze  (t)  =  [e(t)  and  zr  (t)  =  [77(f)  (/)]  .  Note  that  (37)  is  bounded  as  follows 

X\z{t)f  <V{y,t)<X2\z(t)f  (40) 

where  2,,  2,  arc  positive  constants.  After  taking  the  first  time  derivative  of  (37),  the  following  expression  can  be 
obtained 


V  =  -aee2  -  (Ke  +  ae  )se  +  seNe  -  a,rp  -  (Kr  +  a,.  )sr  +  srNr  (4 1 ) 

where  (9),  (18),  (29),  and  the  first  time  derivative  of  (36)  were  utilized.  To  facilitate  the  subsequent  analysis,  the 
following  inequalities  can  be  developed  from  (12)  through  (14)  (see  Appendix  C  for  further  details) 

K<Pl  |ZJ,  K<,p2  \\zr\\  (42) 

where  pt,p2  are  positive  constants,.  By  utilizing  (42),  and  the  triangle  inequality,  V(t)  can  be  upper  bounded  as 
follows 


V<-a$-a\rj( -(Ke+ae)sef  +\se\pl\\ze\\-(Kr +ar]srf  +\sr\p2\\zr\\.  (43) 

By  utilizing  (39),  V(t)  of  (43)  can  be  upper  bounded  as  follows 

V  <  -23||z||2  -Ke\se\2  +  \se\Pl\\ze\\-Kr\sr\2  +  \sr\p2\\z,.\\ .  (44) 

where  23  =  min {ae,ar}.  After  completing  the  squares  for  the  last  four  terms  on  the  right-hand  side  of  (44),  the 
following  expression  can  be  obtained 

y s  -•*,  wr + jjrht + J-wr  )wr  <«> 
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where  A,  =  min  , — —[ .  Provided  that  K,  and  K,.  are  selected  so  that  the  condition  A,  >  A,  is  satisfied,  then 

[4Ke  4Kr\  e  r 

the  following  inequality  can  be  developed 


V  <  W(y )  (46) 

where  W ( y )  e  SJ1  denotes  the  following  non-negative  function 

W(y)  <  -y||z||2  (47) 

where  y  denotes  a  positive  constant.  From  (37),  (40),  (46),  and  (47),  it  is  clear  that  V{y,t)&Lx\  hence, 
e(t),r/(t),se(t),sr(t),ze(t),ze(t),z(t),y(t)GLm.  From  (4),  and  Remark  1,  it  is  clear  that  T;(t),Tr(t)  e  Lx .  From  (9),  it 
is  clear  that  e(t),  ij(t)  e  Lx  .  Thus,  from  Remark  1,  it  is  clear  that  T  (t),  Ty(t)  e  Lr  .  The  previous  bounding  statement 
can  be  used  along  with  (15)  and  (16)  to  prove  that  ue(t),ur(t )  e  Lx .  Thus,  from  (6),  it  is  clear  that  mr(t),mf{t)  e  Lx  ; 
hence,  it  is  clear  that  mc(t),H(t )  e  Lx  .  A  direct  application  of  Theorem  8.4  in  (FI.  Khalil,  2002)  can  be  used  to 
prove  that  z(t)  — >  0  as  t— >oo;  thus,  e(t),se(t\ri(t),sr(t)  — »  0  as  t  — >  oo  ,  where  (39)  was  utilized.  Based  on  the 
definition  of  e{t)  and  r/{t)  in  (4),  it  is  clear  that  if  e{t),rj(t)  — >  0  as  t  — >  oo  ,  then  Te(t)  — >  T;d(t),T  (t)  — >  Trd(t)  as 
t  — >  oo  .  Since  (40),  (46)  and  (47)  hold  globally,  from  Theorem  8.4  in  (FI.  Khalil,  2002),  it  is  clear  that  this  proof 
provides  a  global  asymptotic  stability  (GAS)  result. 

APPENDIX  C 

Upper  Bound  Development  for  Robust  Control 

By  substituting  (13)  and  (14)  into  (12),  the  expressions  of  Nc{te,Te,t )  and  Nr(fr,Tr,t )  can  be  written  as  follows 

Ne  =  Ceaee  +  e,  Nr  =  Cra,:f]  +  rj  (48) 

which  may  be  upper  bounded  as  follows 

Ne  <  Ceae  \e\  +  \e\ ,  Nr  <  Crar  |/;|  + 1^| .  (49) 

Using  the  definition  of  ze(t)  and  zr(t) ,  given  in  (39),  Ne{te,Te,t),  and  Nr(fr,Tr,t )  can  be  bounded  as  follows 

Ne<Pt  Ik  II,  Nr<p2\z,\.  (50) 


Appendix  D 

NOMENCLATURE  LIST 


ae  positive  control  gain 

/?,.  positive  constant  [Rad/sec.m2] 

Ce  engine  block  capacity  [kJ/°K] 

Cpc  coolant  specific  heat  [kJ/kg.°K] 

Cpa  air  specific  heat  [kJ/kg.°K] 

Cr  radiator  capacity  [kJ/°K] 

e  engine  temperature  tracking  error  [°K] 

ea  initial  engine  temperature  tracking  enor  [°K] 


ess  engine  temperature  steady  state  error  [°K] 
s  effectiveness  of  the  radiator  fan  [%] 
i]  radiator  temperature  tracking  error  [°K] 

H  normalized  valve  position  [%] 

mc  pump  coolant  mass  flow  rate  [kg/sec] 
til  f  fan  air  mass  flow  rate  [kg/sec] 
mr  radiator  coolant  mass  flow  rate  [kg/sec] 
M j  pump  coolant  mass  flow  rate  meter 
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M,  radiator  fan  air  mass  flow  rate  meter 
Pt  valve  power  sensor 

P2  water  pump  power  sensor 

P3  radiator  fan  power  sensor 

PM  cooling  system  power  measure  [W] 

P  cooling  system  power  consumption  [W] 

Py  valve  power  consumption  [W] 

pe  positive  constant 

Qin  combustion  process  heat  energy  [kW] 


Qa  radiator  heat  lost  due  to  uncontrollable  air 
flow  [kW] 

sgn  standard  signum  function 

ta  initial  time  [sec] 

7]  coolant  temperature  at  engine  outlet  [°K] 

T2  coolant  temperature  at  radiator  outlet  [°K] 

T}  ambient  temperature  sensor  [°K] 

Te  coolant  temperature  at  the  engine  outlet  [°K] 

T,,  surrounding  ambient  temperature  [°K] 

Tr  radiator  outlet  coolant  temperature  [°K] 

Ted  desired  engine  temperature  trajectory  [°K] 
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